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SUMMARY 



An investigation was made of the effect of the wing 
wake on the lift of the horizontal tail surfaces. In the 
development of expressions for this effect, the growth of 
wing circulation and wing wake, the time interval repre- 
sented by the tail length, and the development of lift 
by the tail were considered. The theory has been applied 
to a specific case to show the magnitude of the effect to 
be expected. 

It is shown that,' for motions below a certain fre- 
quency, the development of . lift by the tail may be repre- 
sented by a simple lag function. The lag is, however, 
somewhat greater than that indicated by the tail length. 



INTRODUCTION 



During unsteady motions of the airplane, the wing 
leaves in its wake a sheet of vortices of varying strength. 
The velocity induced by these vortices may have a pro- 
nounced effect on the direction of the air flow near the 
tail, particularly during motions involving rapid changes 
of lift such as oscillations of short period or passage 
through gusts. 

An approximation to the effect of the wing wake has 
been used by Cowley and Glauert (reference l).« They as- 
sumed that the downwash associated with a change in lift 
is equal to the corresponding steady value but that the 
effect at the tail is delayed by the time required for the 
airplane to travel a distance equal to the tail length. 
It is known, however, that during increases of circulation 
the^wing develops counterrotating vortices which must, for 
a time, at least, induce a strong upwash, increasing the 
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lift of the tail. In addition to the time lag considered 
"by Cowley and Glauert , Doth the .variation of vertical ve- 
locity and the delay in the development of lift "by the tail 
are considered in the present paper. 

The flow around'the wing, hence the wake produced "by 
the wing, is assumed to he uninfluenced hy the presence of 
the tail surface. The interference, is thus confined to 
the effect of the wing on the tail and, since the wake 
formed hy an isolated wing is known, the interference can 
be directly calculated for any relative position of the 
two surfaces. 

Although the theory is thus applicable to a variety 
of arrangements, computations to cover all conditions were 
not considered to be worth while.. In particular, the exact 
vertical location of the tail surface (within the usual 
range) was not expected to be critical. The tail surface 
was therefore considered to be located directly in the 
wake where the effect is a maximum. The effect of tail 
length was investigated and it was found that the results 
obtained from computations covering a typical case could 
be extrapolated to take account of this factor in a sat- 
isfactory manner. 

LIFT FUNCTIONS 



As in other problems involving unsteady flow, it is 
convenient to assume at first a sudden unit jump in the 
angle of attack of the wing and to develop more general 
solutions by operational methods (reference 2). The lift 
of the tail surface under these conditions is due solely 
to • interference from the wing and is to be added to "the 
lift -independently developed. 




>. - 

Relative mo;tion of vortex 



Figure 1. 
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Figure 1 shows the notation -used in the development. 
Variable quantities are expressed in terms of the distance . 

tr t . 

along the flight path, s = -~- , where t represents time. 

c/2 

The analysis is kept nondimen si onal "by expressing all veloc- 
ities in terms of the flight velocity, U 0 , and all lengths 
in terms of the half chord, c/2. 

The rate of development of vorticity by the wing fol- 
lowing a sudden unit increase in angle of attack a 
(a = l(s)) is given in reference . 3 and shows directly the 
distribution of vortices in the wake . The total circula- 
tion at any instant after such a change is denoted by 
r a (s). The vertical velocity of air in the vicinity of 

•the tail induced by a unit wake vortex, Wp(s), (fig. 2). 

may be calculated from the Mot-Savart rule. The result- 
ant variation of vertical velocity follov/ing a unit change 
of angle of attack of the wing follows from the combina- 
tion of these two functions, i.e., ' 

B 

w w (s) = w r (s) r a <0) + J w p (s - s Q ) IV U 0 ) ds Q (1) 

0 

The effect of the vertical velocity, w w , on the tail 
surface may be treated as the effect of a varying gust. 
The lift on an airfoil penetrating this gust is given by 

s 

% (s) = °1 <■> V 0) + / °L (s - s o ) w w'(* 0 ) ds o ( 2 ) 



w 



0 



where C L ( s ) is the gust lift function for the tail sur- 
face (reference 3 ) . 

The integrations (l ) and ( 2 ) may be carried out in a 
single step by using the operational equivalents of the 
functions involved. Thus 

C L, Cb) = ^ (D) l(s) a C L (D ) f a (D) w_(D) l(s) (3) 
l ff t w g J- 

The function 0 L (s) is the lift on the tail surface due 

t w 

solely to a unit change in the angle of attack of the wing 
and is to be added to the lift developed independently by 
the tail surface. 
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The function r( s ) denotes the circulation around 
the wing (or the equal and opposite circulation measured 
around the wake) in the plane of symmetry. The spanwise 
distribution of vorticity is assumed to remain elliptical. 
A unit increment of F(s) thus involves a single wing and 
wake vortex of unit strength at the center, the strength 
falling off toward the tips in accordance with the ellip- 
tic loading. The two vortices are connected by a sheet 
containing only the spanwise component of discontinuity. 
This arrangement can "be derived by the superposition of 
vortices of the type shown in figure 3. 

The centroid of wing circulation is assumed to remain 
stationary at the center of the wing chord. Although the 
wing circulation originates at the trailing edge, little 
error is incurred through this assumption because the trav- 
el of the centroid to the center of the chord is very rapid 
(reference I?). If the wing circulation is replaced by a 
single vortex A (fig. 3), the vertical velocity at the 
tail due to vortices A and B is given by 

w Pa (s) . - 1 - (V~~— - ~-) (4) 
i AB 2tt VI - s l x / 

The operational equivalent of this function is 

D e~ lD E i (ID) - ~] (5) 



AB 



(D) * 



2tt 



where the symbol Ei represents the exponential integral 
function, i.e., 

•~x 



Ei(x) - I --- du 

CO 



The downwash due to the spanwise component of discon- 
tinuity (vortices 0 and D, fig. 3) may be determined 
from figure 6 of reference 3, which shows the downwash at 
the edge of a sheet of discontinuity of varying length. 
The downwash in the region of the tail corresponding to 
any position or extent of, the wake may be obtained by add- 
ing the effects of two sheets of different lengths, as in- 
dicated in figure 4. 
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APPLICATION OP THEORY 



In order to illustrate the application of the theory 
and. to show the order of magnitude of the results to "be 
expected, the lift functions are determined for a specific 
example. The proportions considered are as follows. (See 
fie;. 1.) 

Aspect ratio of wing ----- 6 
Aspect ratio of tail surface -. - 3 

l--„- ________ 5.54 

\ — - 6. 54 

Chord of tail surface - - - - - 1.0 

The effect of the nonci rculatory component of the flow 
about the wing is neglected, its influence at the tail 
being small and constant in value. 



The approximate expressions used for the functions 
C L ( D ) and r a (D) are : 

C L (.33) S 3.77 - _ii56_D_ ^ l_.044_D_ (g) 
* D + .1.116- D + .6.40 

F (D) § 4.71 I^.25_D 0^800_D_ (?) 

a "' D + 0.290 D + 0.690 D + 0,276 

It should be noted that the function C^g is in 

terms of. the half-wing chord which, in the example chosen, 
is twice the corresponding dimension of the tail surface 
itself. There is a slight change in the expression as 
given in reference 3 to make it more closely approximate 
the starting value. 

If the proper values are substituted in equations (4) 
and (5), the function w-p "becomes 

AB 

w ( s ) . ± ( 1 LJ) ( 8 ) 

F AB ^5.54 - s 6.54/ 

and 
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No concise expression for the spanwi se component of 
the downw a sh, w„ , can be given. The curve can he ap- 

1 CD 

proximated, however, by a series of exponential terms of 
the form Ae . For the proportions considered 

w p (s) ~ - 0.083 + 0.047e~ 0 * 067S + 0 . 145e"°- 3 54S 



1 CD 



- 0.156e~°' 742S + 0.047e~ 1 - 45S (10) 



l r (D) S - 0.083 + + ... (11) 

1 CD . D + 0,057 

These expressions are to be added to equations (8 ) and (9), 
respectively, to give the function Wp(D) required in the 

evaluation of equation (3). 

The calculation of C L according to equation (3) 

t w 

results in the lift of the tail surface, as a function of 
the distance traveled following a sudden unit jump in the 
angle of attack of the wing, i.e., the "indicial lift"; 
it is shown in figure 5. 

For a unit change in angle of attack of the airplane 
as a whole, the lift developed independently by the tail 
must be added. This lift increment is given by equation 
(29) of reference 3 for an airfoil of aspect ratio 3 but 
it must be expressed in terms of the wing chord. The 
function C L ( s ) (fig. 6) shows the lift resulting from 

the unit change of angle of attack of the entire airplane. 

Although the indicial lift curves (figs. 5 and 6) 
show infinite values, it is to be noted that the integra- 
tion of the expressions by superposition for any probable 
disturbance results in finite lift at all points. Fur- 
thermore, when the effects of moderate rates of change in 
the angle of attack are integrated, the exact form of the 
indicial response curve is not critical. 
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This point is "best illustrated "by integrating the re- 
sponse to a continuous oscillatory variation of angle of 

attack. If a = e ins , 

C L . (s)=-3 T (D) e ins = -5 T (in) e ins = [a (n )+iB (n ) 1 e ins 
L (a=e ins ) • L J 

The real part, A, is the component of the response in 
phase with the disturhance* and the imaginary part, B, 
is the component 90° out of phase with the disturbance. 

The evaluation of the exponential integral with an 
imaginary argument is given in reference 4, page 80. 



Ei(in) = Ci(n) + i 



[«<») + |] 



where Oi and Si are, respectively, the cosine-integral 
and the sine-integral functions. 

Figure 7 shows the A and the B components Of the 
oscillatory lift function, 0 L (in), the lift on the 

tail surface induced "by a continuous sinusoidal oscilla- 
tion of the wing. Tor a continuous vertical oscillation 
of the airplane as a whole (changes of angle of attack 
without rotation), the function 0^ (in) (fig. 8) shows 
the resulting lift. a 

If the oscillatory lift functions are approximated "by 
a Fourier series, this series will "be found to correspond 
to an approximation of the indicial lift function in the 
form of steps. Thus the function C L (in) is closely 

V 

approximated, as shown "by "broken lines in figure 7, "by the 

expression . 

0 L (in) r 0.30 - 2.20 e -7-i 4 * n (12) 
t w . 

for values of n less than 0.35. If the argument (in) 
is replaced "by the operator D, the resulting, function 
is the operational equivalent of a simple step function:, 
which is an approximation of the corresponding indicial 
lift function. Thus 

0 Lt ( s > 5 0.30 - 2,20 e- 7 » 14D | l(s) £ 

w «- J r i 

0.30 - 2.20 l(s ~ 7.14) (13) 
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This approximation is shown by the broken line in figure 5, 

Similarly, the approximation shown by the broken line 
in figure 8 is 

0 L+ (in) £'3.85 - 1.98 e - 7>48in (14) 
a 

and 

C L (s)~ |3. 85-1, 98 e" 7 ' 485 ] l(s) £ 3.85-1.98 l(s - 7.48) 

'° L J (15T 

and is shown by the broken line in figure 6, 

In the case of an airplane executing pitching motions 
durin? which the angle of attack of the wing does not 
change, the component of response out of phase with the 
disturbance is insignificant. The response may therefore 
be considered instantaneous. Although the general case in* 
volves motions that combine changes of angle of attack, 
a, and of angular displacement, Q , the lift increments 
resulting from each motion, have been separately treated. 
In this form, the results are directly applicable to the 
differential equations of motion of the airplane. 



OONCLUDIKG REMARKS 



Although the indicial lift of the tail surface actu- 
ally shows a 'pronounced variation, it is permissible to 
consider the effect of a simple lag if only moderate rates 
of motion are involved. The lag ,f uncti ons shown in fig- 
ures 5 and 6 give ..the lift quite' accurately during any mo- 
tion that can be compounded of f requenci es lower than one 
cycle in 18 half, chords (n < 0.35). ' ," . ' ' ' , 

It will be noted that these expressions-' differ from 
those assumed by Cowley and Glauert in two ways. Pirst, 
the value of the function from 0 to V is not zero but 
is a positive value, which accounts for the up wash that 
the tail initially . encounters . Second, the distance after 
which the value of the function becomes negative is not 
equal- to, the tail length, I, but occurs at a distance 
somewhat greater than the tail length. This distance ac- 
counts for the lag in the. growth of downwash and the lag 
in the development of lift by the tail and may be called 
the effective tail length. ' - ■ 
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Although the results thus far have "been obtained for 
one specific wing and tail arrangement, the effective tail 
length and the magnitude of the effect to be expected may 
he determined for other cases. 

According to the theory, the pattern of the wake 
formed by the wing remains unchanged as it passes down- 
stream. In addition, the rapid changes of induced veloci- 
ty at the tail occupy a fairly short distance immediately 
ahead of and back of the edge of the wake. For any usual 
tail length, therefore, the time history. of the lift on the 
tail is not substantially altered in relation to the in- 
stant at which the wake strikes the tail.- This point is 
illustrated by figure 9, where the induced vertical veloc- 
ity, wp, following a unit change in wing circulation 
(r = l(s)) is shown for different positions of the wing 
ahead of the tail surface. The principal effect of a 
change in tail length is to shift the origin along the 
relatively flat portion of the indicial lift curves shown 
in figures 5 and' 6. This change in tail length then cor- 
responds to an equal change in effective tail length. The 
effective tail length for any case is thus determined. 

If. the wing wake passes either above or below the 
tail .surface, the lift" functions will not show infinite 
values as they do in the case considered. The peak value 
of the lift function is lowest when the wake passes below 
the tail surface. The final value of the lift is but lit- 
tle affected by the vertical displacement as long as this 
displacement is small relative to the span. 

The theory may be extended to show the effects of 
vertical gusts. The lift of the tail surface due to in- 
terference from the wing during penetration of the gust 
may be calculated with the aid of the curves of figure 5, 
provided that the variation of the angle of attack of the 
wing is known. A close approximation 'to this variation 
of angle of attack during penetration of a varying gust is' 
obtained by measuring the angle of attack with respect to 
the relative wind direction at a point one-fourth of the 
chord ahead of the trailing edge of the wing. This approx- 
imation is based on a well-known result of the thin-airfoil 
theory and is valid as long as the rate of change of gust 
velocity along the flight path is less than that repre- 
sented by an oscillation of one cycle in 18 chord lengths. 
The lift independently developed by the action of the gust 
on the tail surface may be determined from reference 3. 
This lift is directly added to that developed by interfer- 
ence from the wing.. 
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As noted earlier, only in cases of extremely rapid 
changes in wing lift (n > 0.35) is the exact form of the 
indicial lift curve important. Such changes may occur, 
however, in sharp gusts and it is believed that the exten- 
sion of the present investigation to cover the effect of 
vertical position: of the tail surface in these cases would 
"be worth while. 



Langley Memorial Aeronautical Laboratory, 

National Advis-ory Committee for Aeronautics, 
Langley Field, 7a., July 23, 1940, 
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Figure 2.- Variation of vertical velocity at the horizontal tail following a unit 
increase in wing circulation, wj-(s). 
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Figs. 5,6,9. 
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Figure 5.~ Variation of lift on horizontal tail following an Instantaneous 
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Figure 6.~ Variation of lift on horizontal tail following an instantaneous 
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a unit iaorsaM ia wing circulation, tr r (s). 
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Fig. 7. 
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figure 7.- Influence of wing oscillation on the lift of the hori- 
zontal tail. Ct (in) =A(n)+iB(n) 
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Fig, 8. 
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Figure 8.- Lift functions for horizontal tail of airplane in verti- 
cal oscillation without pitching. Cr (in) = A(n) + iB(n) 



